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Regularity of Gaussian Processes on Dirichlet spaces 
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Abstract 

We are interested in the regularity of centered Gaussian processes [Zx{ui))x^m indexed by 
compact metric spaces (M, p). It is shown that the almost everywhere Besov space regularity of 
such a process is (almost) equivalent to the Besov regularity of the covariance A(x, y) = E(ZxZy) 
under the assumption that (i) there is an underlying Dirichlet structure on M which determines 
the Besov space regularity, and (ii) the operator K with kernel K(x,y) and the underlying 
operator A of the Dirichlet structure commute. As an application of this result we establish 
the Besov regularity of Gaussian processes indexed by compact homogeneous spaces and, in 
particular, by the sphere. 

Heat kernel, Gaussian processes, Besov spaces. 
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1 Introduction 

Gaussian processes have been at the heart of probability theory for very long time. There is a huge 
literature about it (see among many others [30], [27], [28] [2], [1] [32]). They also have been playing 
a key role in applications for many years and seem to experience an active revival in the recent 
domains of machine learning (see among others [34], [36]) as well as in Bayesian nonparametric 
statistics (see for instance [45], [24]). 

In many areas it is important to develop regularization procedures or sparse representations. 
Finding adequate regularizations as well as the quantification of the sparsity play an essential role 
in the accuracy of the algorithms in statistical theory as well as in Approximation theory. A way 
to regularize or to improve sparsity which is at the same time genuine and easily explainable is to 
impose regularity conditions. 

The regularity of Gaussian processes has also been for a long time in the essentials of probability 
theory. It goes back to Kolmogorov in the 1930s (see among many others [18], [42], [44] [26], [29]). 

In applications, an important effort has been put on the construction of Gaussian processes 
on manifolds or more general domains, with the two especially challenging examples of spaces of 
matrices and spaces of graphs to contribute to the emerging field of signal processing on graphs 
and extending high-dimensional data analysis to networks and other irregular domains. 

Motivated by these aspects we explore in this paper the regularity of Gaussian processes indexed 
by compact metric domains verifying some conditions in such a way that regularity conditions can 
be identihed. 

In effect, to prove regularity properties, we need a theory of regularity, compatible with the 
classical examples: Lipschitz properties and differentiability. At the same time we want to be able 
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to handle more complicated geometries. For this aspect we borrow the geometrical framework 
developed in [14], [23]. 

Many of the constructions for regularity theorems are based on moments bounds for the in¬ 
crements of the process. Our approach here is quite different, it utilizes the spectral properties 
of the covariance operator. In particular, we use the Littlewood-Paley theory (this point of view 
was implicitly in [12]) to show that the Besov space regularity of the process is (almost) equivalent 
to the Besov regularity of the covariance operator. Especially, it is shown that the almost every¬ 
where Besov space regularity of such a process is (almost) equivalent to the Besov regularity of the 
covariance K{x,y) = K^Z^Zy). 

It is also important to notice that unlike many results in the literature, the regularity is expressed 
using the genuine distance of the domain, not the distance induced by the covariance. 

We illustrate our approach by revisiting the Brownian motion as well as the fractional Brownian 
motion on the interval. We show the standard Besov regularity of these processes but also prove 
that they can be associated to a genuine geometry which finally appears in a nontrivial way. 

We also illustrate our main result on the more refined case of two points homogeneous spaces 
and the special case of the unit sphere in 

In the two subsequent sections, we recall the background informations about Gaussian processes, 
the geometrical framework introduced in [14], [23], and how it allows to develop a smooth functional 
calculus as well as a description of regularity. In Section 4, we state and prove the main result of the 
paper. Section 5 details the case of the standard Brownian motion and fractional Brownian motion. 
In this case, the salient fact is not the regularity result (which is known) but the original geometry 
corresponding to these processes. Section 6 deals with positive and negative definite functions on 
two points homogeneous spaces. Section 7 establishes the Besov regularity of Gaussian processes 
indexed by the sphere. 

2 Gaussian processes: Background 

In this section we recall some basic facts about Gaussian processes and establish useful notation. 

2.1 General setting for Gaussian processes 

Let (n. A, P) be a probability space. A centered Gaussian process on a set M is a family of random 
variables Zx{u}) with x € M and w € 11 such that for all n G N, xi,...,Xn S M, and ai,... ,an G M 

n 

aiZxi is a centered Gaussian random variable. 

i=l 

The covariance function K{x,y) associated to such a process {Zx)xeM is defined by 

K{x,y) :=lK{ZxZy) for {x,y) G M x M. 

It is readily seen that K{x,y) is real-valued, symmetric, and positive definite, i.e. 

K{x,y) = K{y,x) gM, and 

Vn G N, V xi,..., x„ G M, Vai, ..., G M, aiajK{xi, xj) > 0. 

ij<n 
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Clearly, K{x, y) determines the law of all finite dimensional random variables • • •, Z^^)- 

Conversely, if K{x,y) is a real valued, symmetric, and positive definite function on M x M, 
there exists a unique Hilbert space H of functions on M (the associated RKHS), for which K is 
a reproducing kernel, i.e. f{x) = (/, ^(x, •))h, V/ G H, Vx G M (see [5], [37], [15]). Further, if 
is an orthonormal basis for H, then the following representation in H holds: 

= '^Ui{x)ui{y), 'ix,y £ M. 
iei 

Therefore, if is a family of independent A^(0,1) variables, then 

Zxiuj) ■= '^Ui{x)Bi{u}) 
i&I 

is a centered Gaussian process with covariance K{x,y). Thus, this is a version of the previous 
process Zx{u:). 

2.2 Gaussian processes with a zest of topology 

We now consider the following more specific setting. Let M be a compact space and let /r be a 
Radon measure on (M, B) with support M and B being the Borel sigma algebra on M. Assuming 
that (H, A, P) is a probability space we let 

Z : {M,B) (8) (H,yl) i— )■ Z^ioj) G M, be a measurable map 

such that {Zx)x&M is a Gaussian process. In addition, we suppose that K{x,y) is a symmetric, 
continuous, and positive definite function on M x M. Then obviously the operator K defined by 

Kf{x):=[ K{x,y)f{y)dy{y), f£L‘^{M,y), 

Jm 

is a self-adjoint compact positive operator (even trace-class) on L‘^{M,y). Moreover, K{Lp‘) C 
C{M), the Banach space of continuous functions on M. Let > 1^2 > • • • > 0 be the sequence 
of eigenvalues of K repeated according to their multiplicities and let {uk)k>i be the sequence of 
respective normalized eigenfunctions: 

/ K{x,y)uk{y)dy{y) = VkUk{x). 

JM 

The functions Uk are continuous real-valued functions and the sequence {uk)k>i is an orthonormal 
basis for Lp‘{M,yi). By Mercer Theorem we have the following representation: 

K{x,y) = '^i^kUk{x)uk{y), 

k 

where the convergence is uniform. 

Let Ti C L‘^{^,P) be the closed Gaussian space spanned by finite linear combinations of 
{Zx)x£M- Clearly, interpreting the following integral as Bochner integral with value in the Hilbert 
space P, we have 

Bk{uj) = —= / Zxiu})uk{x)dy{x) G P. 

Jm 
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Furthermore, 


and 


^{Bk) = E 


E{BkBi) = 




1 


'M 


Zx{uj)uk{x)d^x{x'^ = / E{Zx)uk{x)d^x{x) = 0 




1 


\/^\/^ Jm Jm 
1 


y/^y/n 


e ( / Zy{uj)uk{y)d^{y) / Z^{uj)ui{x)dyL{x)) 
^Jm Jm ' 

E (Zj; (w) (w)) Ufc (y) (x) (?/) (x) 

K{x,y)uk{x)ui{y)dn{x)dfi{y) = 


M JM 


1 if = / 

0 if / /. 


As the Bk’s belong to the Gaussian space 7i, Bk is a sequence of independent A^(0,1) variables. It 
is easy to see that 


Zy - y^^E{ZyBk)Bk 


LHP) 


= 0 Hy ^ M. 


( 1 ) 


Indeed, clearly E(Z^) = K{y,y) and 


E{ZyBk) =E(Zy^= [ Z^{uj)uk{x)dfi{x)) =-^ [ K{y,x)uk{x)dfi{x) = y/u^Uk{y). 

^ V Afc Jm ' y/Pk Jm 

Hence 

E(Z2) = K{y,y) = Y,Pkul{y) = Y,{E{ZyBk))\ 
k k 

which implies (1). As a consequence, the process 

Zx{l}) ■■= '^y/u^Uk{x)Bk{uj) 
k 

is also a modification of Zx{oj), i.e. P{Zx = Z^,) = 1, Vx G M. 

We are interested in the regularity of the ’’trajectory” x G M i—>• Za;(a;) for almost all a; G H and 
for a suitable modification of Zx{oj). In fact, we will focus on the version Zx{oj). 


3 Regularity spaces on metric spaces with Dirichlet structure 

On a compact metric space (M, p) one has the scale of s-Lipschitz spaces defined by the norm 

ll/llup. := + 0 <S< 1 . ( 2 ) 

x^y p[x,yr 

In Euclidian spaces a function can be much more regular than Lipschitz, for instance differen¬ 
tiable at different order, or belong to some Sobolev space, or even in a more refine way to a Besov 
space. For this purpose, we consider metric measure spaces with Dirichlet structure. This setting 
is rich enough to develop a Littlewood-Paley theory in almost complete analogy with the classical 
case on see [14, 23]. In particular, it allows to develop Besov spaces Bp^ with all set of indices. 
At the same time this framework is sufficiently general to cover a number of interesting cases as 
will be shown in what follows. We next describe the underlying setting in detail. 
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3.1 Metric spaces with Dirichlet structure 

We assume that (M, fi) is a compact, connected measure space, where // is a Radon measure 
with support M. Also, assume that A is a self-adjoint non-negative operator with dense domain 
D{A) C Let Pt = e~^^, t > 0, be the associate self-adjoint semi-group. Furthermore, 

we assume that A determines a local and regular Dirichlet structure, see [14] and for details [19] 
, [40] ,[38] ,[39] ,[41] ,[11], [21], In fact, we assume that Pt is a Markov semi-group {A verihes the 
Beurling-Deny condition); 

0 < / < 1 and f £ imply 0 < Ptf < 1, and also = 1 m (equivalently AIm = 0). 

From this it follows that Pt can be extended as a contraction operator on U’{M, //) for 1 <p < oo, 
i-e- \\Ptf\\p < WfWp, and PtPsf = Pt+sf, t,s>0. 

The next assumption is that there exists a sufficiently rich subspace D C D{A) (see [11]) such 
that f £ D £ D{A). Then we define a bilinear operator “square gradiant” T : D x D 

by 

r(/, g) := -^[A{fg) - fA{g) - gA{f)]. 

Then T{f,f) > 0 and A{f)gdp = g)dp (Integration by part formula). 


Main assumptions: 

1. Let 


p{x,y) := sup f{x)-f{y) fox x,y £ M. 

r(/./)<i 

We assume that p is a metric on M that generates the original topology on M. 


(3) 


2. The doubling property: Denote B{x,r) = {y £ M : p{x,y) < r}. The assumption is that 
there exists a constant d > 0 such that 

{x, 2r)) < 2'^{x, r)), Vx € M, Vr > 0. (4) 

This means that {M,p,p) is a homogeneous space in the sense of Coifman and Weiss [13]. 
Observe that from (4) it follows that 

p,{B{x,\r)) < cq\^ p{B{x,r)) for x £ M, r > 0, and A > l,co = 2'^. (5) 

d is a constant playing the role of a dimension. 

3. Poincare inequality: There exists a constant c > 0 such that for all x £ M, r > 0, and f £ D, 

inf [ if - A)^d/r < cr^ [ T{f,f)dp. 

J B{x,r) 

As a consequence the associated semi-group Pt = e~^^, t > 0, consists of integral operators of 
continuous (heat) kernel pt{x,y) > 0, with the following properties: 
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(a) Gaussian localization: 
ciexp{-^^g^} 


< Pt{x,y) < 


C3exp{-^^^} 


\lp{B{x, Vi))fi{B{y, ^/t)) yjp{B{x, y/t))iJ.{B{y, y/t)) 

(b) Holder continuity: There exists a constant At > 0 such that 


for x,y € M, t > 0. 


\pt{x,y) -Pt{x,y')\ < Cl 


p{y,y')V exp{-^^} 


^ yjp{B{x, Vi))p{B{y, Vt)) 


for X, y,y' £ M and t > 0, whenever p{y, y') < \/i. 

(c) Markov property: 

/ pt{x,y)dp,{y) = 1 for X G M and t > 0. 

Jm 

Above ci,C 2 ,C 3 ,C 4 > 0 are structural constants. 


( 6 ) 


(7) 


( 8 ) 


Remark 1 The setting described above is quite general. This setting covers, in particular, the case 
of compact Riemannian manifolds. It naturally contains the cases of the sphere, interval, ball, and 
simplex with weights. For more details, see [14]- 


Notation. Throughout we will use the notation \E\ := p{E) and 1e will stand for the characteristic 
function of H C M. Also || • ||p = || • \\lp ■= || • \\lp{m,ii)- Positive constants will be denoted 
by c, c', Cl, C, C',. • • and they may vary at every occurrence. The notation a ~ 6 will stand for 
Cl < a/b < C 2 . As usual we will denote by N the set of all natural numbers and No := N U {0}. 

Although general the setting described above entails a structure, which in particular allows to 
develop a complete Littlewood-Paley theory. Next, we describe some basic traits of this framework 
(see [14, 23]). For any t > 0 the operator Pt := e~^^ is a Hilbert-Schmidt operator: 

\\e~^^\\HS-=[ [ \Pt{x,y)\‘^diJ.{x)dfi{y) < oo. (9) 

Jm Jm 

The doubling property (4) implies that M being compact is equivalent to Diam(M) < oo as well 
as to p{M) < oo. It is also equivalent to 

/ pi{B{y,r))~^dp{y) < oo for all r > 0. (10) 

Jm 

From the compactness of M and the fact that A is an essentially self-adjoint non-negative operator 
it follows that the spectrum of A is discrete and of the form: 0 < Ai < A 2 < • • •. Furthermore, the 
respective eigenspaces TLx^ := Ker(A — A^ Id) are hnite dimensional and 


L^{M,p) = ^nx,. 

k>l 


Denoting by P-Hx,. the orthogonal projector onto Px,, the above means that / = Ylk>i ™ 

for all / G L'^{M,fi). In addition, 

Af = '£^kPnxJ, 'Jf^D{A), and PJ =PuxJ, ^f ^ L\ (11) 

k>l fc>l 
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In general, for a function g G L°°(M_|_) the operator g{y/~A) is defined by 


g{VA)f := ^ g{^u)PH,, f, ^f ^ (12) 

k>\ 

The spectral spaces Sa, A > 0, associated with vCd are defined by 

Sa := 0 «A.. 

\/Air<A 


Observe that T,\ G C and hence Sa C for 1 < p < oo. 

From now on we will assume that the eigenvalues {Xk)k>i are enumerated with algebraic mul¬ 
tiplicities taken into account, i.e. if the algebraic multiplicity of A is m then A is repeated m times 
in the sequence 0 < Ai < A 2 < • • • • We let {uk)k>i be respective real orthogonal and normalized 
in eigenfunctions of A, that is, Au^ = X^u^- 

Let Yii{x,y) := '^k{x)uk{y), d > 0, be the kernel of the orthogonal projector onto 

Then as is shown in [14, Lemma 3.19] 

n 5 (x,x) ~ |5(x,(5)|“^ (13) 

Further, if N{5,M) is the covering number of M (or the cardinality of a maximal d—net), then 

dim(Sj_) ^ f \B{x,Vt)\~^dfj,{x) ^ N{Vt,M) t > 0. (14) 

J M 

A key trait of our setting is that it allows to develop a smooth functional calculus. In particular, 
if 5 G C°°(M) is even, then the operator g{ty/~A) defined in (12) is an integral operator with kernel 
g{ty/A){x,y) having this localization: For any cr > 0 there exists a constant Co- > 0 such that 

\g{ty/A){x,y)\ <c„\B{x,t)\~^{l+t~^p{x,y))~‘', yx,y G M. (15) 

Furthermore, g{ty/A){x,y) is Holder continuous. An immediate consequence of (15) is that the 
operator g{ty/A) is bounded on 

|| 5 r(tVA)/||p < c||/||p, yfGLF{M), l<p<oo. (16) 

For more details and proofs, see [14, 23]. 

For discretization (sampling) we will utilize maximal 6-nets. Recall that a set A C M is a 
maximal d-net on M (6 > 0) if p{x, y) > 6 for all x,y G X, x ^ y, and X is maximal with this 
property. It is easily seen that a maximal d-net on M always exists. Of course, if 5 > Diam(M), 
then X will consists of a single point. The following useful assertion is part of Theorem 4.2 in [14]. 

Proposition 1 There exist a constant 7 > 0, depending only on the structural constant of our 
setting, such that for any A > 0 and 6 := y/d there exists a 6-net X obeying 

2"^llfflloo < max l 5 (^)l < Halloo, ygGJ^x- (17) 
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3.2 Regularity spaces 

In the general setting described above, the full scales of Besov and Tribel-Lizorkin spaces are 
available [14, 23]. For the purposes of this study we will utilize mainly Besov spaces. 

The Sobolev spaces Wp = Wp{A), k > I, 1 < p < oo, are standardly defined by 

< := {/ € D{Ai) : \\f\\„f := \\f\\, + \\Aif\\, < rx>}. ( 18 ) 

Consequently, the Besov space Bp^ = Bp^, s > 0, 1 < p, g < oo, is standardly defined by interpola¬ 
tion as in [33] 

B;^:={LP,W^)^^^, e:=s/k, (19) 

where (^LP,Wp)g ^ is the real interpolation space between and Wp, see [14]. 

The following Littlewood-Paley decomposition of functions will play an important role in the 
sequel. Suppose G is real-valued, even, and such that supp<l) C [—2,2], 0 < < 1, and 

‘h(A) = 1 for A G [0,1]. Let T(A) := 4’(A) — <I>(2A). Evidently supp'I'nM+ C [1/2,2]. Set 

To := T and Tj(A) := T(2-^A) for j > 1. (20) 

It is readily seen that To,T G C°°(M), To,T are even, suppTo C [—2,2], suppTj n M+ C 
2'^+^], j > 1, and = 1 for A G M+. Consequently, for any / G LP{M,p), 

1 < p < oo, {L°° := C) one has 

/ = ^T,(^/I)/ in LP. (21) 

j>0 

Note that this decomposition also holds for distributions f ^ S', naturally defined in the setting of 
§3.1, see [23]. 

The following Littlewood-Paley characterization of Besov spaces uses the functions Tj from 
above: Let s > 0 and I <p,q < oo. For a function / G U’{M, p) we have 

/ G Bl^ ^ j > 0, with {e,} G P. (22) 

Furthermore, if / G Bp^^, then 11/11 b|^ ~ We refer the reader to [14, 23] for proofs and 

more details on Besov spaces in the setting from §3.1. 

We next clarify the relationship between and Lips. 

Proposition 2 (a) For any 0 < s < 1 we have Lips C B^ 

(b) Assuming that k > 0 is the constant from (7), then C Lips for 0 < s < k. 

This claim follows readily from the results in [14, 23]. 

Remark 2 In the most interesting case n = 1, Proposition 2 yields Lip s = B^ ^ for 0 < s < 1. 

4 Main result 

We consider a centered Gaussian process {Zx)x&M with covariance function K{x,y) := K{ZxZy) as 
described in § 2.2, indexed by a metric space M with Dirichlet structure just as described in § 3.1. 
We will adhere to the assumptions and notation from § 3.1. 


4.1 Commutation property 

We now make the fundamental assumption that K and A commute in the following sense: 

Definition 1 If K is a bounded operator on a Banach space B [K G ^(B)) and A is an unbounded 
operator with domain D{A) C B, we say that K and A commute if K(D{A)) C D{A) and 

KAf = AKf, yfeD{A). 

Remark 3 Let A be the infinitesimal generator of a contraction semi-group Pt- Then K and A 
commute in the sense of Definition 1 if and only if 


KPt = PtK, Vt > 0. 


We refer the reader to [16], Theorem 6.1.27. 

We now return to the covariance operator K and the underlying self-adjoint non-negative oper¬ 
ator A from our setting. In light of Proposition 3 our assumption that K and A commute implies 
that they have the same eigenspaces. 

Recall that the eigenvalues of A are ordered in a sequence 0 = Ai < A 2 < ..., where the 
eigenvalues are repeated according to their multiplicities, and the respective eigenfunctions {uk)k>i 
are real-valued, orthogonal, and normalized in Hence 

f = Y.^f:Uk)uk, yfeL^{M,^i). (23) 

fe>i 

Let (nk)k>i he the eigenvalues of the covariance operator K. Thus we have 

Auk = AfcUfc and Kuk = i^kUk, fc > 1. (24) 

Remark 4 As a consequence of the commutation property of K and A, the operator AK is defined 
everywhere on L^(M,//) and is closed as K is bounded and A is closed. Therefore, AK is a 
continuous operator from {M, p.) toLp‘{M,p). Clearly, 

KAf = '^{f,Uk)Xk’^kUk for f G L'^ and hence supAfci/fc = \\KA\\^l 2 ) < 00 . 

k>l 

4.2 Main Theorem 

We now come to the main result of this article. 

Theorem 1 Let {Zx)x&M be a centered Gaussian process with covariance function K{x,y) := 
K{ZxZy) indexed by a metric space M with Dirichlet structure induced by a self-adjoint operator A 
such that K and A commute in the sense of Definition 1. Then the following assertions hold: 

(a) If the covariance kernel K{x,y) has the regularity described by 

sup ||iL(a:, •)||_B|^ ^ < 00 for some s > 0, 

x&M 
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then the Gaussian process Zx(uj) has the following regularity: For any 0 < a < | 

e for almost all uj (-8“ i C 5“ oo)- 

(6) Moreover, there exists a unique probability measure Qa on the Borelian sets of ^ such 
that the (canonical) evaluation process: 

Vx G M, dx '■ 00 G B^ i i-)- oo{x) 

is a centered Gaussian process of covariance 

K{x,y) = / 5x{oo)5y{oo)dQa{oo) [Ito-Nisio representation]. 

I Ret 

-^00,00 

(c) Conversely, suppose there exists a > 0 such that Zx{oo) G B^ ^ for almost all oo. Then 

sup \\K{x,»)\\b^^ < oo. 

Remark 5 A key point is that in the above theorem the Besov space smoothness parameter s > 0 
can be arbitrarily large, while 0 < s < 1 in the case when the regularity is characterized in terms of 
Lipschitz spaces. 

For the proof of this theorem we need some preparation. 


4.3 Uniform Besov property of K{x,y) and discretization 

Observe that since the covariance function K{x,y) is a continuous positive definite function on 
M X M, then from (24) it follows that 

K{x,y) = '^VkUk{x)uk{y) and > 0. (25) 

k 

We next represent the Besov norm ol K{x,») in terms of the eigenvalues and eigenfunctions of K 
and A. 


Theorem 2 Let s > 0. Then 


sup ||iF(x, 




~ max ■ 


1 'n ^kulix) 

x&M 


k-.V>T<i 


sup 2 '^‘ 

i>i 


sup 

x£M 


fe:2J-l<v/A^<2J 


Vkul{x)y (26) 


Proof. Note hrst that from (22) it follows that (with 'hj from (20)) 


sup||iF(x,*)||Bj^^ 


~ sup2'^'®sup \\'i>j{y/A)K{x,»)\\oo- 

j>0 X 


But, using (25) we have ['itj{'\/~A)K{x,»)){y) = Yhk^ji''/^)^kUk{x)uk{y) and hence, applying the 
Cauchy-Schwartz inequality it follows that 

sup |(4'j(V^)iF(x,*))(?/)| = sup V4'j(^Afc)i^fcWfc(x). 

x,y X 
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Consequently, 


sup||ir(x,*)||i 3 ^^ ~ sup2^*sup (27) 

* ’ i * 

Clearly, from (20) we have 0 < 'I'j < 1, supp 'I'o C M+ C [0, 2], and supp n M+ C [2-^“^, 2-^+^] for 
j > 1. Therefore, 

supy~]To(v^)t/fc^|(x) < sup y~] Vku\{x) and 

X , X 

^ vT7<2 

sup^Tj(v^)i^fcw|(x) < sup ^ I'kulix), j>l- 
"" k * 2i-l<vT7<2J+l 

These estimates and (27) readily imply that the left-hand side quantity in (26) is dominated by a 
constant multiple of the right-hand side. 

In the other direction, observe that by construction 'I'o(A) = 1 for A G [0,1] and Tj_i(A) -|- 
Tj(A) = 1 for A G [2J-^2J], j > 1. Hence 

sup ^ Vkul{x) <suY>'^^o{V^kykUk{x) and 
\/A7<1 k 

sup ^ Ukulix) < sup"^'ifj-l{^/>^)ukul{x)+ sup'^'ifj{^/^^)ukul{x), j > 1 . 

X . , ,_ ^ 1 X 1 

2-’“l<\/Ar<2^ k k 

These inequalities and (27) imply that the right-hand side in (26) is dominated by a constant 
multiple of the left-hand side. This completes the proof. | 

The following corollary is an indication of how the Besov regularity relates with the “dimension” 
d of the set M, which appears here through the doubling condition (4). 

Corollary 1 Let 7 > d and s = j — d. Then 

Vk = sup ||7i:(x,*)||ij^^ < c. 

X ’ 

Proof. If k'k < c(\/Afc) , then using (13) and (5) we get for any j >1 and x G M 

Vkulix) < c2-T'(j+i) ^ 2 ^^) < Y ul{x) 

fc:2i-l<vT7<2J fc:2i-l<v/A7<2J fe:v/A7<2J 

= c2-^m2.(x,x) < c2-T'^'|H(x,2-^)r^ < c2-^(T'-'^). 

A similar estimate with j = 0 holds for all k such that y/Xk < 1. Then the corollary follows by 
Theorem 2. | 

Remark 6 Observe that 

sup Y nkulix) < c2~^^ ^ Y ^k= Y J t^kUk{x)dfi{x) < c2~^^\M\. 

^ fc: 2 J-l<^/A 7 < 2 J A;: 2 J-l<^/A 7 < 2 J k:2J-^<^%:<2^ ^ 
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We will utilize maximal 5-nets on M along with Proposition 1 for discretization. For any j > 0 
we denote by Xj the maximal 5-net from Proposition 1 with 5 := such that 


1-11 


loo < max|5(0| < Halloo, 


(28) 


The following claim will be instrumental in the proof of Theorem 1. 

Proposition 3 We have 


sup V Ukulix) ~ max V 

k-.^<l k-.y/X^<l 


and for any j >l 


sup V i/fcUfc(x)~max V 
^ k-.20-^<^<2i fc:2i-l<vOT<2l 

with absolute eonstants of equivalence. 

This proposition follows readily from the following 
Lemma 1 Let Xj he the maximal 5-net from above with 5 := 'y2~^, j > 0, and let 

H{x,y):= ^ akUk{x)uk{y), where Ofc > 0. 


Then 

maxFf(^,^)< sup |i?(a:,?/)| < 4maxiL(^,0- 

x,y^M 


Proof. Clearly H(x,y) is a positive dehnite function and hence \H{x,y)\ < x)H{y, y), 

implying 

max \H{f,r])\ =maxH{^,^). (29) 

i&Xj 

Evidently, for any fixed x € M the function H{x,y) G Tj 2 j as a function of y and by (28) 

sup \H{x,y)\ < 2 max |Ff(x, r/)|. 

y(zM 

Now, using that H{x,t]) G S 2 j as a function of x, we again apply (28) to obtain 
sup \H{x,y)\ < 2 sup max |Ff(x, ^)| = 2 max sup \H{x, i])\ 

x,yS.M xdM xdM 

< 4maxmax |Ff(f, ? 7 )| = 4maxFf(f,f). 

v&Xji&Xj' ' i&Xj 

Here for the last equality we used (29). This completes the proof. | 
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4.4 Proof of Parts (a) and (c) of Theorem 1 

(a) Assume sup^-g^f ||Ar(x, •)||s^ ^ < oo. Let {Bk{ijj))k>i be a sequence of independent ^"(0,1) 
variables. Then as alluded in S2.2 


:= '^y/u^Uk{x)Bk{uj) 
k 

is also a version of Zx(uj). Let Tj, j > 0, be the fnnctions from (20) and observe that / £ B^ ^ if 
and only if ||/||bj^_^ ~ Ej>o 2^1l^i(^)/lloo < oo. Clearly, 

(^j(\/A)Z,(u;))(x) = '^^j{y/>^)y/T^Uk{x)Bk{uj). (30) 

k 

For each x £ M this is a Ganssian variable of variance 

< c2"^*- 

k 

Here we used that ^‘j{'/Xk) < Ij the assumption sup^-gj^^ || A(x, •)||b|^ ^ < oo, and Theorem 2. 

For any a > 0 we have 

E(j]2^“||T,-(VI)Z.(a;)||oo) = ^ 2^'“E(||T,(VI)Z.(a;)||oo) 

j 3 

~ ^2^“E( sup |(T,(VI)Z.(o;))(0|) 

< + log(card(A,))i/2. 

3 

Above for the equivalence we used (28) and for the last inequality the following well known lemma 
(called the Pisier lemma, see for instance [43], lemma A.3.1): If Zi ,..., are centered Gaussian 
variables (with arbitrary variances), then 

E( max jZfcj) < c(l + log A^)^/^ max (EjZfcp) 


By (14), we have card(Aj) < c2-^'^. Therefore, if a < |, then 

^ 2 ^" 2 -^*/ 2 (i + log(card(A,))i/2 < < oo. 

3 3 

Consequently, E^ 2^^\\ Tj(\/A)Z,( w) ||oo^ < oo and hence x >->■ Zx{uj) £ 0 < a < s/2, 


n;-a.s. 


(c) Suppose now that uj — a.e.,x e-)- Z^iu) £ B^ a > 0. Then by (30) and (22): 


sup 2-^“ IIE^ j (\/Afc) V^Uk {x)Bk (u}) 


< oo, ui — a.s. 
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By (28) this is equivalent to 


sup 2-^“ max 
j 


k 


< oo, 


u) — a.s. 


(31) 


However, {2^" Y.k is a countable set of Gaussian centered vari¬ 

ables. The Borell-Ibragimov-Sudakov-Tsirelson theorem (see e.g. [27], §7), in particular, as¬ 
serts that if {Gt)t&T is a centered Gaussian process indexed by a countable parameter set T and 
saptgT < oo almost surely, then supjg 7 ’E(Gj) < oo. Consequently, (31) implies 

sup E(2^"^Tj(\/A^)x/I^nfc(0-Sfc) < oo. 

Therefore, there exists a constant C > 0 snch that 

But as before, this yields 

™ < max VTg(v^)z^fcWfc(0 

^exo ^eXo 

fc:VA7<l k 


and, for j > 1, 

max z/fcni(0 < 2max^T2_i(Y^)z/fcU^(^)+ 2max^^'2(y^)j,^^2(^^ < ^2-2^". 

^ ' k,2^-^<V>^<20 " k ^ k 

Here we used that Tj_i(A) -|- 'I'j(A) = 1 for A G [2'^“^, 2-^], implying 'I'^_]^(A) -|- 'I'^(A) > 1/2. 

Finally, applying Proposition 3 we conclude from above that sup^^^i < oo. | 


4.5 Ito-Nisio theorem and construction of the Wiener measure 
4.5.1 Wiener measure on associated to K 

We begin with the following 


Lemma 2 Assume s > 0 and I < p < oo, and let Tj, j > 0, be the functions from (20). Then 

f G Y < oo and ||/||b-^^ ~ Y 

j>o ' ’ j>0 

Proof. From (21) we have for any f € 


j>0 


( 32 ) 


implying 


B, 


p.i 


<E,>oll^i(^)/ll 


B. 


p,i 
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For the estimate in the other direction, note that by (22) 

l>0 

However, supp^j n M+ C [2-5'“^, 2-^+^], j > 1, and hence '^e{''/A)^j{VA) = 0 if \i — j\ > 1. 
Therefore, 

j-i<e<j+i 

On the other hand, by estimate (16) it follows that ||'I'o(vy4)(7|L < cllglL, Vg G L^, and hence 
\\'^'i{^/A)'^j{^/A)f\\p < c\\'^fj{^/A)f\\p, implying 

\\^,iVA)f\\B;^^<c2n\^,{VA)f\\p j;||T,(VI)/b»,, <cJ;2^-^||^I/,(VI)/||p<c||/||b»^^. 

j>0 j>0 

The proof is complete. | 

We now precise Theorem 1, (a) with the following 


Proposition 4 (Ito-Nisio property.) 

Under the hypotheses of Theorem 1 and with the functions 'hj, j > 0, from (20), if •)lls^ oo ^ 

oo, then 

e(^||T,-(VI)Z.(o;))||b.^^) ~E(^2^“||T,(VI)Z.(a;)||oc) <oo, (33) 

i>o j>o 

the map 

I : CO y^^'i)j{-\/A)Z,{Lo){-) G 

j 

is measurable, the serie is normally convergent in and the image probability Q on 

satisfies: 

CO G ^ u:ix) 

is a centered Gaussian process with covariance K{x,y). 


Proof. The equivalence (33) follows by the proof of Theorem 1, (a) and Lemma 2. 

As is well known, for any Banach space B with a measure space (H,H), if G is a finite set 
of indices bi G B and Xi{uj) are real-valued measurable functions, then tv e-;- Aij(w)6i is 

measurable from H to B. Hence, 

a; G H I— )■ 'kj(\/A)Z,(a;) = j{^f^)y/vy.Uk{-)Bk{io) G Bf^^^ 

k 

is measurable. Consequently, by almost everywhere convergence 

/ : w G H 1-^ ^ 'I'j(-\/A)Z,(a;)(-) G Bf^ i 
3 

is also measurable, and I*{P) = Q is a probability measure on the Borelian sigma-algebra, such 
that under Q the family of random variables dx 

^ e B^,i ^ w(x) 

is a centered Gaussian process with covariance K{x,y) = f^a u:{x)uj{y)dQ{u}). | 

We next take on the uniqueness of Q. 
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4.5.2 Gaussian probability on Banach spaces. Proof of Theorem 1 (b) 

For details in this section we refer to [10]. 

Let F' be a Banach space and let B{E) be the sigma-algebra of Borel sets on E. Let E* be its 
topological dual, and assume is a vector space of real-valued functions defined on E, and E) 
is the sigma-algebra generated by F. F = is the vector space of continuous bounded 

functions on E, then ^{Cd{E,M.),E) = B{E) is the Borel sigma-algebra. 

As is well known the sigma-algebra 'y{E*,E) generated by E* is B{E) if E is separable (By 
separation B{E)) is generated by open balls and by separation and Hahn-Banach theorem open 
balls are in 'y{E*,E)). 

Proposition 5 Let E be a separable Banach space. Let H be a vector subspace of E*, endowed 
with the a{E*,E) topology. Then 

H is closed H is stable by simple limit. 

Proof. The implication ^ is obvious. We now prove As is a separable Banach space, this 
is a consequence of: 

Vi? > 0,B(0,i?) = {f £E* : ||/||e* < R} is metrizable for a{E*,E) 

and, by Banach-Krein-Smulian theorem, H is £^)-closed if and only if Vi? > 0, B{0,R) D H 

is Fi)-closed. But this is clear, since we only have to verify that for every sequence (fn) C 

B{0, R) Cl H such that lim„_j.oo fn{x) = f{x), Vx G E, we have / G B{0, R) n H, which is what is 
assumed. | 

- cr(E* ,E) 

Corollary 2 If E is a separable Banach space and H is a vector subspace of E*, then H 
coincides with the smallest vector space of functions on E, stable by simple limit containing H. 
Moreover, 

Proof. Clearly, as E* is stable by simple limit (by Banach-Steinhauss theorem), the smallest vector 
space of functions on E, stable by simple limit containing H is contained in E*. And if j{H,E) 
the sigma-algebra generated by H, the vector subspace of E*, 7 (Lf, Li)—measurable is also stable 
by simple limit. | 


Lemma 3 Let E he a separable Banach space. Let H be a subspace of E* separating E. Then 


j{H,E)=-f{E*,E)=B{E). 


There is at most one probability measure P on the Borel sets of E such that, under P, 7 G Lf is a 
centered Gaussian variable with a given covariance K{'y,'y') on H. 

Moreover if such a probability exists, then 


1 . 


E* is a Gaussian space, and E 


:LRe,p) 


is the Gaussian space generated by H. 


2. There exists a > 0 such that 



e"ll^ll£(iP(x) < 00. 


(34) 
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Proof. By the Hahn-Banach theorem = E* and 

-i{H,E)=^{E\E)=B{E). 

Now, if K{'y,'y') is a positive dehnite function on E[, it determines an additive function on the 
algebra of cylindrical sets related to H: 

{x & E : (71 (x),... , 7 n(x)) G C}, ji & El, C Borelian set of M"'. 

Now, the sigma-algebra generated by this algebra is the Borelian of E. 

Assume that such a probability P exists . Let EL = E* D Clearly 

Gaussian space generated by H, and if ( 7 n)n>i G ^ is such that Vx G E, lim„i_^oo 7 n(a^) = l{x) 

_ p\ 

exists, then clearly 7 G by the Banach-Stheinhauss theorem, and -y £ H ’ since a simple 
limit of random variables in a closed Gaussian space belongs to this Gaussian space. Therefore, 
•y £ EL, which by Proposition 5 implies that EL is closed. But H G EL and H ’ = E* leads to 

EL = E*. 

Finally, (34) is just the Fernique theorem. | 


Corollary 3 Let M he a set and let E be a separable Banach space of real-valued functions on M. 
Let 


yt£M, f £E% f{t) £ M. 


If 6t £ E*, then 

B{E) = y{n,E), n = {Y,ai^u}. 

finite 

Let K{x,y) be a definite positive function on M x M. There is at most one probability measure 
P on the Borelian sets of E such that, under P, {6x)xeM is a Gaussian process, with covariance 
K{x,y), and E* is a Gaussian space. 


We now come to the main assertion here. 


Theorem 3 (Wiener measure) In the setting defined above, if K(x,y) a continuous positive 
definite function on M such that sup^-g^^ ||A'(x, ^ < 00 and the associated kernel operator 

K commutes with A, then there is a unique probability measure Q on the Borelian sets of B^.^, 
a < |, such that the family of random variables: 

Vx G M, cj G ^ w(x) G M 
is a centered Gaussian process of covariance K{x,y). 

This theorem holds due to the previous result and the fact that B^ ^ is separable. It also proves 
Part (b) of Theorem 1 . 
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4.6 Regularity and reproducing kernel Hilbert spaces 

Let K{x,y), {x,y) G M x M, be a continuous real-valued positive definite kernel on a compact 
space M. It is well known that K determines a real Hilbert space Mx of functions, for which the 
evaluation: 

Vx e M, 5x ■ f G !->• f{x) is continuous. 

Moreover, 

y K{x,y) = K^{y) G Mk, V/ G Uk, S^{f) = {f,K^)ujt, and is a total set inHi^. 

The space Mx is the completion of span{iL(x, •) : x G M}, more precisely 

■= {Hv) = '^oiiK{xi,y) : ||/i||ii = ^ aiajK{xi,Xj) = 

ieF i,j&F jeF 

It is well known (see e.g. [15]) that 

||h||^ = 0 for h G e?, ^ h{y) = 0, Vy G M. 

It is also well known (see [31]) that 

K{x,y) = '^giix)gi{y) gi G Mx, Vi and {gi)iei is a tight frame for Mx- 

i&I 

In our geometric framework, where the regularity spaces are linked to a suitable symmetric 
positive operator, and K{x,y) is compatible with the geometry, we have 

K{x,y) = '^UkUk{x)uk{y) = '^[^/i^Uk{x)][y^Uk{y)], 
k k 

where {uk)k>i is an orthonormal basis for Lp‘{M,g) consisting of eigenfunctions of A associated to 
the eigenvalues (Afc)fc>i. Therefore, clearly (yT^Ufc)fceN,i/fc^o is a tight frame of H. 

Moreover {5x)x&m C is dense in in the weak cj(]H[)^, Hi^) topology. In fact, the following 
theorem holds. 

Theorem 4 (a) Let N(z^) := {k G N, 7 ^ 0} and define 

= {/ : M^M : f{x) = E Uk{x), (ak) G with inner product 

k&n(v) 

{f,g)H = (^ afcV^Mfc(-), ^ := ((ofe), (/3fc))£2(N(,.))- 

Then TL is a Hilbert space of continuous functions and {^/FkUk)keN{|y) orthonormal basis for 

TL . In fact Mx = TL . 

(b) We have for s > 0 

Mx ^ sup ||iL(x,«)||Bs^^ < 00 . 

xGM 
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(c) Let a < I and denote by J : Mk W = ^ the natural injection, and J* : W* i—)• 

Then 

/m(J(Hx) = spanjufc € VF : fc e N(z^)}. 

Under the probability Qa on W = W* is a Gaussian space, and 

[ ( 7 ~iV( 0 ,||J*( 7 )llk)- 

Jw ^ 

Moreover, W*^ isometrically isomorphie to ~ '^k- 

Proof, (a) We have 


\ak\V^\uk{x)\ < (^ ^ (l«fcP)"( X] ^k\uk{x)\^Y = \\a\\p(j^^^,,))yjK{x,x). 
km{u) fceN(i/) fceN(i/) 

Therefore, "H is a space of continuous function, ^ is a Hilbert space, and {y/nkUk)k&i[y) is an 
orthonormal basis for TL. Furthermore, 

Kx{y) = K{x,y) = ^ VkUk{x)uk{y) = ^ [y/^Uk{x)][y/T%Uk{y)] G TL 

feeN(y) fceN(i/) 


and 

{f,Kx)H = f{x), \\Kx\\n= {y/iTUk{x)f = K{x,x). 

k&n{v) 


Also, clearly, % = EIx- 

(b) Suppose that sup,j-g^ ||Ar(x,< oo and let f{x) = EfceNH where 

(ofc) G f'. Then 

Ttj{'/A)f{x) = ^ ^j{^/^)ak^/iTkUkix), 
km(u) 

implying, for j > 1, 


Tj(VA)/(x)| < ( \ak\YY \'^jiV^)\‘^^k\ukix)\Y^ 

fceN(i/) fceN(y) 

<II/IIhk( Y ^k\uk{x)\Y^ 

k-.2i-^<\k<2i+^ 


where for the last inequality we used the assumption and Theorem 2. Similarly o{y/A)f {x)\ < 

c||/||h^- Therefore, in light of (22), 


<c||/||hk- 


( 35 ) 


Suppose now that (35) holds. Then for every sequence (a^) G with ||(afc )||£2 < 1 we have 
I Y Uk{x) <c2~^^^‘^, Vx G M, 

keN(u) 
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which by duality implies 


^ \uk{x)\^j ^ < c2 3 > 0. 

km{y) 

Just as in the proof of Theorem 2 we get for j > 1 

^ Vkul{x) < ^ ^ |«fc(a;)P < c2 

fc:2J-l<vT^<2J fceN(!^) fceN(i/) 

and similarly Y.k-.^<i^kul{x) < c. Consequently, sup^-g^ \\K{x,»)\\b‘^^^ < oo. 


-JS 


(c) Clearly = spanju^ ^ W : k ^ N(i^)} and from the previous results, under the 

probability Qa on VC = VC* under Qa is a Gaussian space. Let now F C M he finite and 

^ construction 


'W 


^ai6xiico)j dQa{uj) = ^ aiajK{xi,Xj) = 
ieF i,j&F ieF 


We obtain the result by density of the span of {Sx)xeM in I 

Remark 7 Let f G Clearly 

f{uj) : a; G VC = ^ [ f{x)u){x)dfi{x) 

JM 

belongs to VC*. Hence, under Qa,f is a Gaussian variable and 

iE(/)^ = [ ( [ f{x)uj{x)dn{x)) dQa{uj) = [ [ f{x)uj{x)dfi{x) [ f{y)uj{y)d^i{y)dQa{ov) 
^JM ^ JwJm Jm 


iw 



M JM 


f{x)f{y)(^J^(x{x)u;{y)dQa{u})^d^i{x)dfi{y) = {Kf, f) 


Consequently, 


and ||J*(/)||^|^ = {K f, f) 

5 The examples of Brownian Motion and fractional Brownian mo¬ 
tion 


Here we illustrate our main result (Theorem 1) on the example of the standard Brownian motion. 


5.1 Wiener representation of Brownian motion 

Assuming M = [0,1], consider the kernel 

K{x,y) = X Ay = ]^{x + y -\x - y\). 
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It is easy to find the eigenfunctions and eigenvalues of the operator K with kernel K{x, y). Indeed, 
we wish to find sufficiently smooth solutions cj) of the problem 


X A yc/){y)dy = \(l){x), where (/>(0) =0, A / 0. 


Differentiating both sides of the above identity we obtain implying = 0. An¬ 

other differentiation leads to (j)''{x) + X~^4> = 0. As a result, we obtain the following eigenfunctions 
and eigenvalues: 

(j)k{x) = V2 sin {k + l^)7rx, Xk = ^ ^ , /c = 0,1,- 


Therefore, 


Kix,y) = 


The associated Gaussian process takes the form 


sin 


(7r(A:-b 2 ))^ 


[k + hTTxsin [k + 7:)T^y- 


ZxW) = Yj 71-fT— A/2sin {k -|- -Jttx ■ Bk(u;), where Bk ~ A^(0,1), i.i.d. 

A natural Dirichlet space (with Neumann-Dirichlet boundary conditions) is induced by the 
operator 

Af:=-r, D{A) :={/€C2]0,l[nCM0,l] : /(O) =/'(I) = O}. 

Clearly, 

/ Af{x)f{x)dx = [ f'‘^{x)dx, f e D{A), 

Jo Jo 


and 


1 1 2 1 

A(sin(A;-|- -)vr«)(x) = {{k + -)7r) sin (/c-|- -)7rx. 


Also, the distance on [0,1] is defined by 

p{x,y)= sup \f{x)- f{y)\ = \x-y\. 

In this setting, the Poincare inequality and the doubling property are obvious, and clearly 
\K{x,y) - K{x,y')\ <\y - y'\, implying sup ||A:(3 :, •)||bi^ ^ < 1. 

xSM 


(36) 


So far everything looks fine, unfortunately the Dirichlet space induced here does not verifies all 
the conditions described §3.1, e.g. the associated semi-group is not Markovian due to the fact that 
the function 1 ;= 1m does not belong to D{A). 

In the next subsections we will discuss a useful way to circumvent this problem, in particular, 
we will identify a Dirichlet space adapted to the framework of Brownian motion. This will require 
careful study of positive and negative definite kernels. 
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5.2 Positive and negative definite functions 

For this subsection we refer the reader to [7], [35], [9]. Recall first the definitions of positive and 
negative definite functions: 

Definition 2 Given a set M, a real-valued funetion K{x, y) defined on M xM is said to be positive 
definite (P.D.), if 

n 

K{x,y) = K{y,x), and Vai,..., q;„ G M, Vxi,..., G M, aiajK{xi,Xj) > 0. 

*,i=i 

As shown in §2.1 the following characterization is valid: 

A:(x,y) is P.D. K{x,y) =¥.{Z^Zy), 

where {Zx)xeM is a Gaussian process. 

For any u € M we associate to K(x,y) the following P.D. kernel 

Kuix,y) := K{x,y) -i-K{u,u) - K{x,u) - K{y,u) = E[{Zx - Zu){Zy - Z„)]. 

Clearly, 

Ku = K K{u,u) = 0. 

Definition 3 Given a set M, a real-valued funetion 'ip{x, y) defined on MxM is said to be negative 
definite (N.D.), if 

'ip{x,y) ='ip{y,x),\/x,y & M, fi{x,x) = 0, and 

n 

Vai,...,On G M s.t. Aj = 0, Vxi ,... ,Xn G M, aiajfi{xi,Xj) < 0. 

i i,j=l 

The following characterization is valid (see e.g. [7, Proposition 3.2]): 

ijj{x,y) is N.D. fi{x,y) =E{Zx- Zyf, 

where {Zx)x&M is a Gaussian process. Consequently, ip{x,y) > 0, Vx,y G M. 

From above it readily follows that ^/'tjj{x, y) verifies the triangular inequality: 

|v/'0(x,y) - yjfi{z,y)\ < 'dx,y,zGM. (37) 

The following proposition can easily be verified. 

Proposition 6 (o) Let K{x,y) be a P.D. kernel on a set M, and set 

'4’K{x,y) := K{x,x) + K{y,y) - 2K{x,y). (38) 

Then fix is negative definite. The kernel fix will be termed the N.D. kernel associated to K. In 
faet, if K{x,y) = E{ZxZy), then fix{x,y) = E{Zx — Zyf^. Furthermore, fix = fixuj ^ ^ ■ 

(b) Let fi be a N.D. kernel, and for any u £ M define 

N{u, fi){x, y) := ]^[fi{x, u) + fi{y, u) - fi{x, y)]. 

Thus, if fi{x,y) = E{Zx — Zy'fi, then N{u,fi){x,y) := ^[(^ 3 , — Zu){Zy — Z„)]. Then N{u,fi) is 
P.D. Moreover, 

N{u,fix) = Ku- 
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The next assertion contains our key idea. 


Proposition 7 Let 'ip{x,y) he a real-valued continuous N.D. function on the compact space M, 
and set 


K{x,y) 


1 

2|M| 



u) + u) - ipix, y)]dy{u). 


Then 

(a) K is positive definite. 

(b) 1 is an eigenfunction of the operator K with kernel K{x,y), that is, 


/ K{x,y)l{y)dy{y)= K{x,y)dy{y) = Cl, with C = —^ 
Im Jm 


il;{u,y)dp{u)dp{y) > 0 . 


M JM 


(c) 


3z G M s.t. K{z,z)=0 K{x,y) = 0 


fi{x,y) = 0 . 


Proof. Parts (a) and (b) are straightforward. For the proof of (c) we first observe the obvious 
implications: 

'fi{x, y) = 0 => K{x, y) = 0 => K{z, z) = 0, Vz G M. 

Now, let K{z,z) = 0 for some z G M. Then 2 |^ ~ ^)]^/^(^) = 0- 

By definition 'il}{z,z) = 0 and hence Jj^fi{z,u)dfi{u) = 0. However, fi{z,u) is continuous and 
fiiyZ, u) > 0. Therefore, 'il>{z, u) = 0, Vu G M. Now, employing (37) we obtain for x,y G M 

V'iP{x,y) = \^/'^p{x,y) - y/'ijj{z,y)\ < y/fi{x,z) = 0, 

and hence 'fi{x,y) = 0. This completes the proof. | 

Remark 8 The following useful assertions can he found in e.g. [7], [35], [9]. 

Let fj{x,y), defined on M x M, obey fj{x,y) = fi{y,x), Vx,y G M, and ip{x,x) = 0. Then 

-0 is N.D. Vt > 0, e-^^ is P.D. 

^[> is N.D. =► V 0 < a < 1, is N.D. 
is N.D. log(l + fi) is N.D. 


The following proposition can easily be verified. 


Proposition 8 Let M be a compact space, equipped with a Radon measure y. Assume that K{x,y) 
is a continuous P.D. kernel and let ip := 'fix be the associated to K{x,y) N.D. kernel, i.e. ip{x,y) := 
K{x,x) + K{y,y) -2K{x,y). Set 

Ku{x,y) ■.= ]^[i){x,u) +fi{y,u) -fi{x,y)\, 

and 

•= FuTT / +T/>(y,n) - V'(x,y)](iy(n) = — / Ku{x,y)dy{u). 

I Jm I Jm 

Denote by K and K the operators with kernels K{x,y) and K{x,y). 
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Then 


(39) 


Kix,y) = K{x,y) + \M\-^ Tt{K) - \M\-^Kl{x) - \M\-^Kl{y). 

Moreover, = ip, = K^, and 

kl = Cl, C = Tr{K)--^[ [ K{x,y)dy{x)dy{y) =[ [ ip{u,y)dy{u)dy{y) >Q. 

In addition, 

K = k + const. Kl = C'\ (40) 

and const. = — |M|“^(Tr(iir) — 2 C"). 

Proof. From the respective definitions, we infer 

'ipix, u) + ipiy, u) — ip{x, y) = {K{x, x) + K{u, u) — 2K{x, u)] 

+ [K{y, y) + K{u, u) -2K{y, w)] - [K{x,x) + K{y,y) - 2K{x,y)] 

= 2[K{u, u) - K{x, u) - K{y, u) + K{x, y)] 

and hence 

k{x,y) = K{x,y) + [ [K{u,u) - K{x,u) - K{y,u)]dy{u) 

pw I Jm 

= K{x,y) + |^(Tr(/s:) - Kl{x) - Kl{y)). 

The remaining is a consequence of Proposition 7. | 

Remark 9 Observe that if K{x,y), ip{x,y), and K(x,y) are as in Proposition 8, then 

3z G M, k{z, 2) = 0 k{x, y) = 0 ^^{x, y) = 0 K{x, y) = const. 

Indeed, clearly we have only to show the implication ip{x,y) = 0 K{x,y) = const. However, 
ip{x, y) = 0 ^ K{x,x) + K{y, y) = 2K{x, y) < 2yJK{x,x)^K{y,y), 
implying {^/K{x, x) — \/K{y, y))^ < 0, which leads to K{x,x) = const, and 

K{x,y) = ]^{K{x,x) + K{y,y)) = const. 

Remark 10 Assume that we are in the geometrical setting described in ^3.1, associated to an 
operator A. Just as in ^40, suppose K{x,y) is a P.D. kernel such that the associate operator K 
commutes with A. From ( 4 I) we have Al = 0. 

Moreover, it is easy to see that 

AIm = 0 and dimKer( 2 l) = 1. (41) 

Indeed, the Markov property (8) yields AIm = 0. To show that dimKer(2l) = 1, assume that 
= 0; / S D{A). Then T{f,f) = 0. Assume that f ^ constant. Then f{x) / /(y) for 
some x,y G M, x ^ y. For T{f,f) = 0 we have T{af,af) = 0 for each a > 0. Then by (3) 
p{x,y) > a\f{x) — /(/)|, Va > 0, implying p{x,y) = 00 , which is a contradiction because M is 
connected (see [14])- Therefore, Af = 0 implies f = const, and hence dimKer(2l) = 1. 

Hence 

AKl = KAl = 0. 

However, as dimKer(A) = 1, necessarily K1 = Cl. Therefore, K = k A const. 
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5.3 Back to Brownian motion 

Assume again that M = [0,1] and K{x,y) := x A y = ^{x + y — \x — y\). We will adhere to the 
notation introduced in §5.2. 

The associated to K{x,y) N.D. kernel ip{x,y) = '^Kix,y) (see (38)) take the form 

V>(x, y) = K{x, x) + K{y, y) - 2K{x, y) 
and the induced P.D. kernel Ku{x,y) becomes 


Ku{x, y) = -[^{x, u) + i/jiy, u) - 'ipix, y)] 

= K{x, y) — K{x, u) — K{y, u) + K{u, u) 

= ]^[\x-u\ + \y-u\-\x-y\], uG[0, 1]. 


Thus we arrive at the following P.D. kernel 


K{x,y) := J Ku{x,y)du = -[x^ + (1 - xf + y"^ + (1 - vf - 2|x - y\]. 
Denoting by K the operator with kernel K{x,y), we have 


Kl{x) = i 


1 />! pi pi 

\y — u\dudy + / \x — u\du — / |x— 
0 JO Jo Jo 





Further, using that cos kirydy = 0 for A: G N we get 


iF(cos k'K*){x) = 


1 


1 /-i 


0 JO 


\f\--yv 


\y — u\ cos kny dudy — - I \x — y\ cos kny dy 

/o 


Integrating by parts we obtain 


Sin kiry ^ 

/ \x - y\cosk 7 Tydy = \x - y\ -;— - ( - l[o,d(y) + l[x,i](2/)) 

Jo y='^ Jo 


sin kiry 
irk 


dy 


— ^ J sin kny dy — J sin kiry dy^ 


(vrfc)^ 


\y — u\ cos kny dudy = 

10 JO Jo 

Putting the above together we infer 


2 cos kiTx 1 + (—1)^ 

+ 


By the same token 
"l /■! 



{nkY 


2 cos kiru 1 + (—1) 
(Tik)'^ (vr/c)^ 


du = 


l + (-l)^ 

[irkY 


~ COS kjTX ~ 

K(cosk7r»)(x) = — —VA: G N, and A1 = 1/6. 
[Tiky 
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Observe also that the functions {1} U (\/2cos /c'7r3:)fcgN form an orthonormal basis for L^{0, 1). 

Let H^{0, 1) be the space of the functions / € -L^(0,1) that are two times weakly differentiable 
and /', f" G -L^(0,1). Consider the operator 

Af := D{A) := {/ G H\[0, 1]) : /'(O) = /'(I) = 0}. 


Clearly, 


/ {Af)gdx = [ f'g'dx = [ fAgdx 
Jo Jo Jo 


10 JO Jo 

and hence A is positive and symmetric. In fact, A generates a Dirichlet space, and also 

coskiTX ^ D{A) and ^(cos A:7r«)(x) = (vrA:)^ cos fevrx. A: > 1. 

Let H^[0, 1] := {/ G L^(0,1) ; \f'{u)\‘^du < oo}. This defines a Dirichlet form: 

A,D{A) = {feH\0,l) : |^'/'(xy(x)dx| <c||(^|| 2 , 'icj, e 

[ f'{x)(p'{x)dx = [ Af{x)4>{x)dx 
Jo Jo 


Thus 


and the distance is defined by 


p{x,y)= sup (/>(x) - (/>(y) = |x - ?/|. 

The Poincare inequality is well known to be true in this case. So we are now in the setting 
presented above. K{x, •) is uniformly Lip 1. Therefore, the centered Gaussian process associated 
to K is almost surely Lip a, a < ^. 

The process Yx{uj) = Zx{u}) — Zq{uj) has the same regularity, and ^(|a;| + \y\ — \x — y\) is the 
associated kernel. This is the Brownian motion, with the above regularity. 


5.4 Brownian motion and fractional Brownian motion, through analysis on the 
circle 

The Laplacian on the torus (for instance on M/2Z) is a typical example of an operator generating 
a Dirichlet space with all properties that are required for defining a regularity structure. If we 
represent 1R/2Z by the arc length parametrisation we have a Dirichlet space associated to: 

Af := -/", D{A) := {/ G - l,l[nCM-l,l] : /(-I) = /(I),/'(-I) = /'(I)}, 

j Af{x)g{x)dx = j f'{x)g'{x)dx, 

and the distance is defined by 


Pix,y) = inf {/(x) - f{y) : \f'\ < !,/(-!) = /(I),/'(-I) = /'(!)} = |a: - y| A (2 - |x - y\). 


Clearly, the eigenfunctions of A are (cos A:7ra;)fcgp}Q and {smk7rx)keN- 
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5.4.1 Brownian motion on the circle 


Using the Fourier series expansion, we have, 


Hence, 


I I 1 4 cos(2n + IIttx , , 

I^I = 2-^E (2n + l)2 • -el-l.l]. 

nSN '' ^ 


, , , , , , ,,14 cos(2n + l)7r(x — u) 

p{x, y) = x-y A{2 - x-y) = ---Y, 2 ^ 

2 vr^ (2n + 


nSN 


implying 


1 4 

^ - k - y| A (2 - |a; - y|) = — 2 ^ 


cos(2n + l)7r(x — y) 


TT^ ' (2re + 1)^ 

nSN ^ ' 


TT 


E 


cos(2n + l)7rx cos(2n + l)vr?/ 4 sin(2n + l)7ra: sin(2n + l)7ry 


(2n + l)2 


+ -E 

TT 


(2n + 1)2 


nEN ^ tiGN 

From this it follows that the kernel K{x,y) is P.D. and 

y) = K{x, x) + K{y, y) - 2K{x, y) = \x - y\ A {2 - \x - y\) = p{x, y). 


Therefore, the Gaussian process Zx{ijj)x&[-i,i] associated to K{x,y) is a Brownian field with 
respect to p. Its regularity is Lip a, a < ^ with respect to the metric p. We can now restrict to 
[0,1], p{x,y) = \x — y\ Vx,?/ G [0,1]. Thus, considering Wx := Zx — Z^, restricted to x G [0,1], we 
get the classical Brownian motion with 

Wo := 0, E(IU^ - Wyf = 2|x - y\ 
and we again obtain its regularity as a byproduct. 


5.4.2 Fractional Brownian motion on the circle and on [0,1] 

Let 0 < a < 1. The Fourier series expansion of |x|“ on [—1,1] takes the form 


X = 


0 + 1 


+ 2 > cos kirx / 

h Jo 


cos kiTudu. 


Integrating by parts we get 


f 


cos kTTU du = — 


+ si 

kTT Jq 


sin knu du = — 


a 


Hence, 


0 + 1 


- |xr = 2o^7fc 
fc>i 


{'KkY+^ Jo 

cos kirx 
'' (7r/c)“+i ’ 


rkn 

I ' 


^sinudu. 
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where 


Ik 



sin u du 



1 


sin u du 



sin rt I hit 



sinu 
{u + 


k-l 

du =: ^^{ — lyuj. 
j=0 


Here oq > oi > • • • > 0 and limj^oo = 0- Hence, 7 := limj^oo Ij exists and 0 < 7 < 7 r"+^/(a+1). 
Therefore, 

Kaix,y) = - i\x-y\ A (2 - |x -y|)" = 2a ^ 7 ^ 

E cos kiTx cos kiry + sin fcvrx sin kny 
(7rA:)“+i 

fc>i ^ ^ 


is a P.D. kernel compatible with the Dirichlet structure dehned by the Laplacian on the circle M/2Z. 
Moreover, as 0 < a < 1 


Ka{x, y) - Kaix, y')\ < \p{x, y)" - p{x, y')“| < p{y, y'T- 


Consequently, the associated Gaussian process {Zx{oj))x(^y_i^^ with covariance function Ka{x,y) is 
Lip /3 for /3 < ^ on [—1,1]. If we restrict this process to a; G [0,1] as p{x, y) = \x — y\, Vx, y G [0,1], 
we get a Gaussian process on [0,1] with covariance — \x — y|", and such that 

¥.{Zx-Zyf = \x-yy. 

Hence, the process {Z^ — ■^o)a;e[o,i] has covariance |x|“ + \y\°‘ — \x — y|“ and regularity Lip^d, j3 < ^. 
This is the standard fractional Brownian function. 


Remark 11 If a > 1, then as above 7 ^ = X]o<i<A;-i(“l)'^ + ^ sinndu and hence 7 ^ > 0 

if k is even and 'jk < 0 if k is odd. From this one can deduce that p{x,y)°‘ is not a definite negative 
function on the circle. 


6 Positive and negative definite functions on compact homoge¬ 
neous spaces 

Here we present some basic facts about positive and negative dehnite kernels in the general setting 
of compact two point homogeneous spaces. Then, in the next section, we utilize these results and 
our main Theorem 1 to establish the Besov regularity of Gaussian processes indexed by the sphere. 

6.1 Group acting on a space 

Let (M, p) be a compact space equipped with a positive Radon measure p. Assume that there exists 
a group G acting transitively on {M,p), that is, there exists a map {g,x)&GxMi-^g-x&M 
such that 
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1. h- {g- x) = {hg) ■ x, Mg, h^G, 

2. 3e G G s.t. e ■ X = X, Mx G M (e is the neutral element in G), 

3. Mx,y G M, 3g G G s.t. g ■ x = y (transitivity), 

lMiM{ 9 )f){x)dy{x) = f{g-^ ■ x)d^i{x) = f{x)dg{x) MgGG ,Mf G L^, 

where {'yig)f){x) := f{g~^ ■ x). Hence, {'y{g))g£G is a group of isometry of L^. 

Definition 4 A continuous real-valued kernel K(x,y) on M x M is said to he G-invariant if 

K {9 ■ x,g-y) = K{x,y), MgGG,Mx,yGM. 

If K is the operator on iP' with kernel K{x,y), then K is called G—invariant if ^{g)K = K^{g), 
Mg G G, that is, 

[ K{g-^ ■ x,y)f{y)dfi{y) = [ K{x,y)f{g~^-y)dg{y), Mf G . 

JM JM 

Remark 12 (o) Assume that K{x,y) is a continuous G-invariant kernel, then 

(i) K{x, x) = K{g ■ X, g ■ x) and hence K{x, x) = \M\~^ Tt{K), and 

(ii) 

/ K{x,y)dg{y)= K{x, g ■ y)dg{y) = K{g-^ ■ x,y)dg{y), MgGG, 

JM JM JM 

and hence 1 := 1m is an eigenfunction of K, that is, K{x,y)l{y)dii{y) = C'l(x). 

(6) Suppose K{x,y) is a continuous positive G-invariant kernel, then 

i^K{x, y) := K{x, x) + K{y, y) - 2K{x, y) = 2{C - K{x, y)) = 2{\M\-^ Ti{K) - K{x, y)), 

clearly 'ijj{x,y) is G-invariant and by (40), K{x,y) = K{x,y) + |M|“^(Tr(K) — 2C'). 

(c) Suppose ijj{x,y) is a G—invariant N.D. kernel and consider the associated P.D. kernel K, 
defined as in (39). Then K{x,y) is G-invariant, and 

If ~ 1 

^ uTT / u)dg{u) = Cq and K{x, y) = Cq - -ifix, y). 

\ Jm ^ 

Thus, in this framework there is one-to-one correspondence up to a constant between invariant P.D. 
and N.D. kernels. 

6.2 Composition of operators 

Let K{x,y) and H{x,y) be two continuous kernels on M x M as above, and let K and H be the 
associate operators. The operator K o H \s also a kernel operator with kernel K o H(x,y): 

KoH{x,y)= / K{x,u)H{u,y)dia{u). 

Jm 

Observe that: 


29 



1. If K{x,y) = K{y,x),H{x,y) = H{y,x) then 

KoH{x,y)= / K{x,u)H{u,y)dy,{u) = / H{y,u)K{u,x)dfi{u) = H o K{y,x). 

Jm Jm 

2. If K{x,y) and H{x,y) are G—invariant, then so is iiT o H. Indeed, 

K o H{g ■ x,g -y) = / K{g ■ x,u)H{u,g ■ y)dg{u) = / K{g ■ x, g ■ u)H{g ■ u, g ■ y)dy{u) 


IM 


IM 


/ K{x,u)H{u,y)dij{u) = K o H{x,y). 

Jm 


6.3 Group action and metric 

Assume that we are in the setting of a Dirichlet space defined through a non-negative self-adjoint 
operator on just as in §3.1. Suppose now that, 


i{a)A = ^ 7 ( 5 ), V 5 e G 


or equivalently 

-i{g)Pt = Ptlig), Vt > 0, V 5 G G, 

i.e. Vt > 0, is G-invariant. Clearly r(/i,/ 2 ) is also G-invariant: r(/i,/ 2 ) ='r{^{g)fi,'y{g)f 2 ) 

and the associate metric p{x, y) is G-invariant: 

p {9 ■ x,g-y) = p{x,y), 'JgeG. 

Definition 5 In the current framework, (M, p, A, p, G) is said to be a two point homogeneous spaee 

Vx, y, x', y' e M s.t. p{x, y) = p{x', y'), 3^ G G s.t. g ■ x = x', g ■ y = y'. 

In particular, V(x, y) £ M x M, 3g £ G s.t. g ■ x = y, g ■ y = x. 

Theorem 5 Let (M, p, A, p, G) be a eompaet two point homogeneous space. Then we have: 

1. Any G-invariant eontinuous kernel K{x,y) is symmetrie. 

2. If K{x,y) and H{x,y) are two G-invariant continuous kernels, then K o H = H o K. 

In particular, if K{x,y) is a G-invariant continuous kernel, then KA = AK. 

3. Any G-invariant real-valued eontinuous kernel K{x,y) depends only on the distance p{x,y), 
that is, there exist a eontinuous funetion A: : M 1 —)• M, such that 

K{x, y) = k{p{x, y)), Vx, y £ M. 


This theorem is a straightforward consequence of the observations from §6.2 and the definition 
of two point homogeneous spaces. 

Let now M be a compact Riemannian manifold and assume that A := —Am is the Laplacian 
on M, p is the Riemannian metric, and p is the Riemannian measure. Also, assume that there 
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exists a compact Lie group G of isometry on M such that {M, p,G) is a compact two 

point homogeneous space. For the link with Gaussian processes see: [ 6 ], [20]. 

Let 0 < Ai < A 2 < • • • be the spectrum of —Am- Then the eigenspaces := Ker(AM + Afc Id) 
are finite dimensional and 

L^{M,p) = ^nx,. 

k>l 

Let Phxi^ y) be the kernel of the orthogonal projector onto Then if K{x, y) is a G-invariant 
positive definite kernel we have the following decomposition of K{x, y), which follows from Bochner- 
Godement theorem ([17], [22]): 

K{x,y) = '^T^kPHxSx,y), z^fe> 0 . 
k>0 

7 Brownian motion on the sphere 

In this section we apply our main result (Theorem 1) to a Gaussian process parametrized by the 
unit sphere in This is a Riemannian manifold and a compact two point homogeneous 

space. More explicitly, 

G = SO{d + 1), H = SO{d), G/H = 

The geodesic distance p on is given by 

pii.p) = arccos(C,r/), 

where is the inner product of ^,77 £ Clearly, 

V^, p £ S'^, Vfir £ G, p{g ■^,g-p) = p{^, p), and V^, p £ 3g e G s.t. g ■ ^ = p. 

Thus G acts isometrically and transitively on Furthermore, 

'7^) P} p' £ s.t. p{^, p) = p{^', p'), there exists g ^ G s.t. g ■ ^ and g ■ p = p'. 

Therefore, is a compact two point homogeneous space. 

Let —Agd be the (positive) Laplace-Beltrami operator on As is well known the eigenspaces 
of — Agd are the spaces of spherical harmonics, defined by 

"Haj, := Ker(Agd + Afc/d), \k ■-= k{k + d - 1) = k{k + 21 /), k>0 

One has L^(S'^) = 0fc>o^Afe and the kernel of the orthogonal projector Phx^ onto 77 Aj, is given by 

Pn,^{(,ri) = LUx) := IS^I-l(l + 

Here G^{x), k >0, are the Gegenbauer polynomials defined on [—1,1] by the generating function 

^^^ = y r^Gfc^(x). 

(1 — 2xr + 

^ ’ fc >0 
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Therefore 


E k{k + 2i>)P'n^^f 


k>0 


and the invariant continuous positive definite functions on are of the form 

K{C,r]) =J2’^kLi{{^,T])) = J2’^kLi{cosp{^,rj)), 


where 


Note that 




< oo. 


l: 


;( 1 )|S" 1 = [ L^((C, 0 )rfMO=dim(i^A,(§")) 


k + d 
d 




Let 


nix) := ^ Clearly. W^il) = sup = 1, 


Then (see 


lirnSW 

-0 C£(l) 

li.. 


= Tk{x) (= PT^(x) by convention), 
= (= W^ix) by convention). 


Here is the Chebyshev polynomial of first kind (Tfc(cos0) = cos k6). The invariant continuous 
positive dehnite functions on S'^ are of the form 


K^{C,v) = Y^alWi:i{C,rj)) = Y,^kWk{k^osp{C,v)), 4 > 0, j;4<oo. 


k>0 


k>0 


Clearly, 


'^a’'kWk{cosp{^,p)) = '^-^L’^kicosp{^,p)), L^(1)~A:'^ ^ 

z/fc = 


Therefore, 


(42) 


= O 


dim(^Aj \k^-^J' 

The following Schoenberg-Bingham result (see e.g. [8]) plays a key role here: If f is 
a continuous function defined on [—1,1], then f{{f,,p)) is a positive definite function on and 
invariant with respect to SO{d + 1) for a/1 d G N if and only if 

fix) = Y1 ttnX^, where an > 0 and E“» = /w< 


oo. 


n>0 


n>0 


Therefore, for such a function / 

fix) = J2^kWkix), al>0, and ^ ^ a*, =/(I), 


k>0 


k>0 k>0 


and hence 


fii^,v)) = '^a'j)Wlf{{^,p)) = '^j^^L'j){{C,p)) = /(cosp(^,ry)). 


k>0 


k>0 
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7.1 Fractional Brownian process on the sphere 
Theorem 6 For any 0 < a < 1 the function 

is negative definite, and the associated Gaussian process has almost everywhere regularity 
7< f. 

Proof. Consider first the case when a = 1 (Brownian process). We will show that for some 
constant C > 0 the function C — rf) is an invariant positive definite function. To this end, by 
Schoenberg-Bingham result we have to prove that there exists a function 

f{x) = ^ anX^, with a„ > 0 , '^On < oo, 

n>0 


such that /{{(,, rj)) = f {cos p{f,r])) = C — p{(,,p). Luckily the function | — arccosa: does the job. 
Indeed, it is easy to see that 


fix) 


TT 

— — arccos x = arcsin x = 


E 

j>o 


2h 




( 2 ^ ^2j+l 


and 


E 

j>o 


ilhihh 

jKDj 


^ (Gauss). 


Here we use the standard notation {a)j := a{a + 1) • • • (a + j — 1) = r(a + j)/r(a). Therefore, 

TT TT 

/((?: 7 )) = - - arccos(^, 77 ) = --pif,p). 

Clearly, |/((C,f/)) — fi{f,jh'))\ E pip^p') and by Theorem 1 the associated Gaussian process 
(Z|(w))^g§d is almost surely in (hence in Lips) for 0 < s < ^. Furthermore, 

¥.{Zl - Z^)^ = 2/(1) - 2f{{C,p)) = 2p{C,p). 

Consider now the general case: 0 < a < 1 {Fractional Brownian process). From above it follows 
that ijj{f,,p) ■= p{C,p) is an invariant negative dehnite kernel. Then the general theory of negative 
dehnite kernels yields that for any 0 < a < 1 the kernel p) = p(?) p)‘^ is invariant and negative 
dehnite. Therefore, for a sufficiently large constant C > 0, 

K{C,p) = C-^p{C,p)^ 

is an invariant positive dehnite kernel. On the other hand, 

\K{^,p)-K{C,p')\ = i|p(e,r7)“ -p(C,r7')“| < ^p{p’,p))^. 

By Theorem 1 it follows that the associated Gaussian process (Z|(a;))gggd is almost surely in B^^ 

7 < ^, and hence in Lips, s < ^, and the proof is complete. | 


Remark 13 From the definition of the process, we have 

E(Z|‘-Z“)2 = p(^,7?)“. 

This directly connects to the regularity proofs of such a process using generalization of Kolmogorov- 
Csensov inequalities. See for instance [3] and [25]. 
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7.2 Regularity of Gaussian processes on the sphere: General result 
Theorem 7 Let 

^ a > 0. 


n>0 


ni 




Then 

K{^,rj) := f{cos{^,r])), d>l, 

is an invariant positive definite function, and the associated Gaussian process (u;))^g§d is almost 
surely in ^ for 7 < a. 

Proof. By Corollary 1, it suffices to show that f{x) can be represented in the following form (see 
(42)): 


fix) = Y,BjW^ix), 0<Bj = implying uj = 0 (^ 7 ^) = O(v^) 


\2a-\-d 


By lemma 1 in [8] and the obvious identity r(x + n) = (a;)nr(x) we obtain the representation 

^ n\ 


X = 




0<2fc<n 


Substituting this in the definition of f{x) we obtain 

/w = E#-” = E# E {j = n-2k) 


n>0 


i>o 


n—2k=j 

1 A 


2^k!(n)n-k+i 




j>o 


2 . 22^k!(n),+k+i 


-.^BjW^ix), 


j>o 


where 


Br-= 


{j + n){2n)j 


A 


j+2k 


j\2i 

^^22^fc!(re),+,+i 

if + 

^j+2fc 

j!2i(re)j+i 

^ 2‘^^k\{v + j + l)fc 

{2v)j 

^j+2fc 


^ 22fcA:!(zy + j + l)fc ■ 


However, for n > a we have 

Cl (a) r(n -a) 02 ( 0 ) 


n 


l + Q — 


re! re^+“ 


and 


hence ^ ^ = 0(r(re - a)). 
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We use this to obtain for j > a (with c = c{a)) 


E- 


A 


j+2k 


r(j + 2k — a) 

22fc/c!(i/ +j + l)fc “ 2‘^’^k\{u + j + l)fc 


s^E 


= cr(j -a)^ 


r(j +2k — a) 


k>0 


T{j-a) 2‘^>^k\{i2 + j + l)k 


= cT{j -a)Y^ 


(j - a) 2k 


k>0 


22^ kl{u + j + l)k 


= cT{j - a)'^(^ 
k>0 


j / j - a + l 

2 )k\ 2 


kk\{v + j + l)k 


where we used the Legendre duplication formula (see e.g. [4]): 

(^) 2 fc _ r (6 + 2 A:) _ / 6 \ /6 + 1 \ 

22fc “ 22fcr(6) “ V2 Av 2 )k 


By the Gaussian identity (see e.g. [4, Theorem 2.2.2]) 

n-a + l \ 1 T{i 2 + j + l)r {,2 + j + 1 - 2^ - 2^) 

2 )k\ 2 )kk\{u + j + i)k r(i. + j + i-^)r(i. + j + i-i^) 

_ r(z/+ j + i)r(z/+ i + a) 
r(i/ + i +1 + f )r(z/ + i + i 

5 ^n/I M c'C' 

B < „ (2^)j r(j - a)r(z/ +j + l)r(r/+^ + g) 

- ''j\2i{u)j T{i2 + i + l + q)r{u + i + 1 + f) • 

Applying again the Legendre duplication formula, we get 


We use this above to obtain for j > 2a 

^ ^ {2iy)j r(j - a)T{i' + j + l)r(z^ + ^ + a) 

^ - ''j\{u)j T{\)T{2u + 3 + 1 + a)2-2^-“ 

_ r(2zy + j)r(i/) T{j - a)T{u + j+ l)T{v+ \ +a) 

~ ""+{3 + l)r{ 2 u)T{i 2 + j) r(i)r(2i/ + j + l + a)2-2^-“ 

= c2"+^f? + + ^ + r(j-a) _ T{2u + j) 

’ r{u + ^) T{j-a + l + a)T{ 2 i 2 + j+ l + a) 

1 _ 1 c 

— ^ 3 + k> ^ (+^l+a (^2i/ + j)l+“ “ jl+2a ■ 


Here we used once again the the Legendre duplication formula. It is easy to show that Bj < c(a), 
if 3 < 2a. Therefore, Bj = 0( ^i+2a ) and this completes the proof. | 
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Corollary 4 Let a > 0, b > 0, c> a-\-b, a = c — a — b, and let 


n 


{a)n{h)n 
{c)n n! ■ 


Then Fafi-dieh)) invariant positive definite function on the sphere S'^ and the associated 

Gaussian process has regularity 7 < a, almost everywhere. 
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